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In a collinear antiferromagnet with easy-axis anisotropy, symmetry dictates that the spin 
wave modes must be doubly degenerate. Theses two modes, distinguished by their opposite 
polarization and available only in antiferromagnets, give rise to a novel degree of freedom to 
encode and process information. We show that the spin wave polarization can be manipulated 
by an electric field induced Dzyaloshinskii-Moriya interaction and magnetic anisotropy. We 
propose a prototype spin wave field-effect transistor which realizes a gate-tunable magnonic 
analog of the Faraday effect, and demonstrate its application in THz signal modulation. 

Our findings open up the exciting possibility of digital data processing utilizing antiferro¬ 
magnetic spin waves and enable the direct projection of optical computing concepts onto the 
mesoscopic scale. 

Spin waves are propagating spin precessions in magnetically ordered media. Since spin waves 
can carry pure spin currents in the absence of electron flow, they are considered to be potential 
information carriers for low-dissipation, spin-based computing technologies, known as magnon- 
ics [1-4]. The possibility of using waves instead of particles for computing also enables new device 
concepts for data processing, such as spin wave logic gates [5, 6]. As a first step towards magnon- 
ics, it is necessary to encode binary data into spin waves. Similar to electromagnetic waves, spin 
waves are characterized by their amplitude, frequency, and polarization. However, in ferromagnets 
(FMs), the spin wave polarization is always right-handed with respect to the background magneti¬ 
zation. Therefore, one usually has to adopt the spin wave amplitude [7] or its frequency to digitize 
information. 

By contrast, it is well known that the spin wave dynamics of a collinear easy-axis antiferromagnet 
(AFM) admits two degenerate modes with opposite circular polarization [8, 9], as illustrated in 
Figs. 1(a) and (b). These two modes can be recombined into an equivalently orthogonal but linearly 
polarized basis, as shown in Figs. 1(c) and (d). This two-fold degeneracy places AFM spin waves 
in a similar situation as electromagnetic waves, which suggests that the spin wave polarization can 
be harnessed to encode information. In fact, since the magnon chirality is connected to the photon 
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polarization, optical methods have been exploited to excite the two degenerate modes [ 10 , 11 ]. 
According to recent investigations [12-14], selective excitation and detection of the two circularly- 
polarized modes are also achievable via electron spin currents with corresponding polarizations, 
paving the way to encode information into the polarization of AFM spin waves. 

The next crucial step towards AFM magnonics is to control this degree of freedom via external 
fields in order to perform logical operations on the encoded data. Since the degeneracy is protected 
by the combined symmetry of time-reversal (T) and sublattice exchange (X), a viable control must 
resort to interactions that break either or both of the two symmetries, i.e., an effective field that 
couples to the Pauli matrices <t = {ui, cj 2 , < 73 } spanning the doubly generate space. 

In this article we propose that the Dzyaloshinskii-Moriya interaction (DMI) [15-17] can be used 
for such a purpose. The DMI is expressed generically as T>ab ■ {Sa x Sb), where T>ab is the DMI 
vector that couples two spins Sa and Sb] in an AFM, Sa and Sb represent the two antiparallel 
moments in a magnetic unit cell. Since the DMI changes sign upon sublattice exchange {A -H- B), 
it breaks the degeneracy between the two circular modes. We show that the DMI in AFMs behaves 
as an fictitious field that couples to in the degenerate space, leading to opposite phase shifts for 
the two circular modes. As a result, when a linearly-polarized spin wave is subject to the DMI, 
the opposing phase shifts of its circular components lead to a rotation of the linear polarization 
direction, which realizes a magnonic analog of the Faraday rotation of electromagnetic waves. If 
we identify the x- and y-polarized spin wave modes as 0 and 1 in binary operations, a rotation by 
7 r /2 then corresponds to a NOT operation in magnonic computing. 

Based on the Faraday rotation of AFM spin waves, we propose a gate-tunable field-effect tran¬ 
sistor serving as the magnonic analog of the Datta-Das device of electrons [18]. We demonstrate 
its application in the amplitude-shift keying as a THz signal modulation. Finally, by including the 
field-induced anisotropy [19], we can realize direct transitions between the two circularly-polarized 
modes, which, together with Faraday’s rotation controlled by DMI, enable a complete control 
of spin wave states over the entire Bloch sphere. Possible electrical detections of the spin wave 
state on the Bloch sphere are also discussed. Our findings open up the possibility of digital data 
processing harnessing antiferromagnetic spin waves, and enable the direct projection of optical 
computing concepts onto the mesoscopic scale. 

Results 

Spin wave spectrum. Under the continuum description, a collinear AFM is characterized by the 
staggered field n = {Sa — Sb)/‘2S and the small magnetization m = {Sa + Sb)/‘^S. Consider a 
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quasi-one dimensional nanostrip with an easy-axis along the 2 -direction, as schematically shown 
in Fig. 2(a). A perpendicular gate voltage is applied to break the inversion symmetry and induce 
a DMI of the following form (see Supplementary SI) 


Hd = 2?[n • (V X n) -|- V • (n X m) — m • (V x m)], 


( 1 ) 


where T) is the DMI strength and V = y x V with y being the mirror plane normal. The DMI 
may be nonzero even without gating, as the geometry of the interface already breaks the mirror 
symmetry. In the exchange limit, m <C n, we will drop the last term; the second term is a total 
derivative that does not affect the local dynamics. This leaves us with only the first term of Eq. (1). 
The total action in terms of n(x, t) is [20] 


S= dtdx 




dn 


dt 


, r- 2 , duy dn^ 


(2) 


where A = Ja/2 is the stiffness with J > 0 the antiferromagnetic exchange coupling and a the 
lattice constant, c = 2aJ/h is the spin wave velocity, and K. = K/i2a) represents the easy axis 
anisotropy. For small DMI, the (classical) ground state of Eq. (2) is the uniform Neel state. If the 
DMI exceeds a threshold value Dth = 2\/}CA = VKJ, the ground state twists into a spiral. But 
for our purposes, all discussions are restricted to the sub-threshold regime where the uniform Neel 
ground state is preserved. 

Linearize Eq. (2) in terms of the small deviation n± = {nx,ny} of the staggered field from its 
equilibrium value = z, and define '4>± = nxdziuy to associate with the right-handed (left-handed) 
mode. Setting 6S/5ijj± = 0, we obtain a two-component Klein-Gordon equation 


= [{AV^ - /C) + cJ3D(-fV)]T , 


where 'k = (^/;+,^_)^ . Solving the above equation with the ansatz ip± = 'ip±e 




( 3 ) 


yields the 


dispersion 


wl = c^[k^ ^Qk + Z], 


( 4 ) 


where Q = T) j A and Z = K/A. Equation (4) is plotted in Fig. 2(b). For an arbitrary given 
frequency a; above the antiferromagnetic resonance point at ojr = c\/Z, we find that the splitting 
of wave vectors 


A/c = k+ — k- = Q (5) 

is independent of to. Let L be the length of the gated region, then the linear polarization of an 
AFM spin wave will rotate by 7r/2 when AkL = vr, which can be regarded as the magnonic analogy 
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of the Faraday rotation of an electromagnetic wave. This mechanism also has a direct analogy 
in the electron spin field-effect transistor (FET) [18], where AA: is frequency independent and is 
proportional to the Rashba spin-orbit coupling. 


AFM spin wave field-effect transistor. The key to realize an AFM spin wave field-effect 
transistor is the gate-tunable DMI [17, 21]. It was shown that in transition metal compounds, an 
electric field E generates an DMI of the strength T) = JaeE/Ego, where £so originates from the 
spin-orbit interaction and is typically on the order of 3 eV [21]. The electric field E can be produced 
by either a gate voltage or the interfacial mirror symmetry breaking. Recently, a gate-tunable DMI 
has been observed in a ferromagnetic insulator [22]; but the same mechanism is applicable to AFMs 
as well. The condition AkL = tt mentioned above is tantamount to an electric field 


^ Ego 
~ 2 ~^' 


( 6 ) 


If L ~ I/xm, Eq. (6) is satisfied by E ~ 3 V//im. The required E field can be scaled down by 
increasing L, but care should be taken that L not exceed the magnon phase coherence length 
ifj,. Though the actual value of depends on multiple factors, it suffices to consider the magnon 
dephasing due to the Gilbert damping. This allows for a simple estimate ~ Vp/{aoj), where Vp 
is the phase velocity, a is the Gilbert damping constant, and w is the spin wave frequency. Since 
Vp = oj/k, we have £(j, ~ l/{ak). Eor a ~ 0.005 and k < 2/im“^, we find that exceeds 100 /rm. 
This indicates that spin waves in AFMs can propagate over large distance without losing phase 
information and provides a loose upper limit for L. On the other hand, the need to suppress spiral 
formation sets a lower limit for L — maintaining the stability of the Neel ground state yields a 
maximum allowed electric field 

( 7 ) 

ea \ J 

which depends on the ratio K/ J that can differ by several orders of magnitude in different mate¬ 
rials (0.1 ~ 10“^) [23-26]. The value of E^ax then sets the corresponding lower bound on L via 
Eq. (6), which can be as short as a few nanometers. 


THz data modulation. Since the resonance frequency of AFMs is typically in the THz range, 
the AFM spin wave FET has important applications in high frequency data modulations. As an 
example, we perform a numerical simulation of the amplitude shift keying (ASK) based on the 
device schematic of Fig. 2(a). Assume that the spin wave is generated at one end of the chain 
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by an oscillating magnetic field along the y-direction at 1.4 THz. The spin wave is subsequently 
modulated by a 20 GHz square wave through the gate voltage Vg, which covers 330 atomic sites 
(see Method). Fig. 2(c) plots the y-component of the transmitted spin wave, where we see a clear 
on/off ratio. 

While the bit rate of such an ASK is as high as THz, the actual signal transfer rate (Baud rate) 
is limited by the relaxation time r of the spin wave FET. As indicated by the dotted red line in 
Fig. 2(b), T is the transient period that a transmitted spin wave adapts to an abrupt change of 
the gate voltage, so the maximum Baud rate cannot exceed 1/r. Physically, the relaxation time 
can be regarded as the time required for the wavefront of an incident wave to traverse the gated 
region: r must therefore scale as Ljc. Since c = 2aJlfi, we have r ~ {L/a)uj~^^ with ojj = J/h. 
In typical AFMs, ojj is around hundreds of THz, so r is only a few picoseconds when L/a ^ 100. 
This allows reliable signal transfer at a Baud rate below 0.1 THz. The capacity of signal transfer 
can in principle be enlarged by reducing L as long as the ground state remains collinear. 

Figure 2(d) shows the spatial pattern of the staggered field. Long after the relaxation time, 
the trajectories of the Neel order on every 60 atomic sites (30 unit cells) are plotted for several 
periods. While the transmitted wave is essentially linearly-polarized, the reflected wave and the 
wave inside the gated region slightly open up into elliptical precessions. This fact implies that 
the relative ratio between the left-handed wave and the right-handed wave is close to 1 in the 
transmitted region, whereas it slightly deviates from 1 in the reflected and gated regions. By a 
straightforward wave matching calculation at the gate boundaries, we find that the deviation of 
this relative ratio from 1 is linear in T)j J (= 0.5% in our simulation) in the reflected and gated 
regions, while it is proportional to ( 2 ?/J)^ in the transmitted region, which explains the observed 
pattern (see Supplementary S2). 


Spin wave state on the Bloch sphere. Manipulations of AFM spin waves are not limited to 
the Faraday rotation, which is realized by the coupling of us via the DMI. Coupling to other Pauli 
matrices extends the attainability of spin wave states to the entire Bloch sphere (see Fig. 3). It has 
been shown that besides the DMI, an electric field could also induce magnetic anisotropy if it is 
tilted towards the easy axis [19]. This amounts to applying an artificial magnetic field along e\ in 
the internal degenerate space, which rotates the spin wave state in the e 2 — 63 plane, toggling the 
spin wave chirality. In contrast, the DMI studied above behaves as a magnetic field along 63 that 
rotates the spin wave state along the equator. Under the basis of the circular modes, the equation 
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of motion for spin waves is given by {A/(P‘)d‘f'^ = with the effective Hamiltonian 


H = 





cost'd Q + a3'Dsm'd^{—iV), 


( 8 ) 


where is the polar angle of the gate voltage direction with respect to the jz-axis, and K,' is the 
strength of the i?-field induced anisotropy. This Hamiltonian enables us to explore the entire Bloch 
sphere spanned by the wave function [cos 0 / 2 , sin 0 / 2 ]"'^, where 0 and (j) are the spherical angles 
specifying the position on the Bloch sphere as shown in Fig. 3. 

The mapping of the AFM spin wave onto a Bloch sphere places it on an equal footing as the 
electron spin, in the sense that it can store information through coherent superpositions of two 
states. Therefore, the AFM spin wave can be regarded as a classical analog of the quantum bit. 
However, unlike the challenging task of manipulating quantum information, manipulating AFM 
spin waves is much simpler. For example, by tunning the direction of the gate voltage (i.e., the 
angle '0^ in Eq. ( 8 )), the artificial fields that couple to ai and <73 could be made equal, which 
mimics the behavior of a Hadamard gate in quantum computing. 


Detection of spin wave state. Besides conventional optical approach [10, 11], a spin wave 
state {9,4>) can also be read off by virtue of spin pumping [14]. Suppose that a heavy metal 
is deposited in direct contact to the AFM at the far end of a spin wave FET [see Fig. 2 (a)]. 
Spin waves transmitted from the gate can pump spin current into the heavy metal, and this spin 
current is converted into the inverse spin Hall voltage [1, 27] that is monitored by two voltmeters 
Vi and V 2 [28]. By measuring the DC component V 2 and the effective AC component Vi (the 
root-mean-square value), one can determine 0 and (j) by (see Supplementary S3) 


D2 

= arccos:^, 


(f> = 2 arcsin 


Vi 

\/ 2 E™ V 2 




(9a) 


where is the maximum DC voltage along x (for purely right-handed spin wave), and 17™ is the 
value of Vi at zero gate voltage (also its maximum). In principle, the effective AC component can 
be determined by measuring the output power of the circuit. 


Discussion 

In contrast to FMs, AFMs are devoid of the long range dipolar interaction. In thin film FMs, 
the dipolar interaction substantially modifies the spectrum in the long wavelength limit [29]. In 
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particular, the wave vector shift AA: due to the DMI tends to zero near the FM resonance fre¬ 
quency [30, 31]. In AFMs, however, AA: is a constant proportional to the DMI strength as long as 
a; is greater than the resonance freqnency, as shown above by Eq. (5). Therefore, the DMI induces 
the change of wave vectors more efficiently in AFMs than FMs. 

Our proposal for the AFM spin wave FET hinges on the ability to generate sufficiently strong 
DMI via a gate voltage, which can be realized in compounds with heavy elements. Another crucial 
requirement is the availability of suitable AEMs with only easy-axis anisotropy. Possible candidates 
include MnF 2 [23], FeF 2 [24], and RbMnEs [25]; all are antiferromagnetic insulators. However, the 
Neel temperatures of these materials are too low (below 100 K) for room temperature applications. 
A promising replacement is NaOsOa with the Neel temperature around 410K [26]. It is a G-type 
collinear AEM with a dominating easy-axis. To fully unlock the potential of AFM spin wave FETs 
thus calls for further development of room-temperature easy-axis AFMs. 

Methods 

The simulations were performed by solving the Landau-Lifshitz-Gilbert equation on a ID chain 
using the Dormand-Prince method. We set the saturated staggered magnetization to be unity and 
scaled everything with frequency. The parameters were then taken as ujj = 100 THz, = 500 
GHz, and uok = 10 GHz. In order to avoid effects on the gate system due to reflection from the 
chain boundaries, we made the length of the chain much larger than the gate system. 

The full chain was constructed with 10® atomic sites (5 x 10® AFM unit cells). At the center of 
the chain, a local magnetic field of 1.4THz was used as a source of linearly polarized spin waves. 
Near this source, a 330 site region was gated with DMI. To verify our predictions, we observed the 
spin wave profiles exiting the far side of the gate. 
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Right-handed mode 

(a) 



FIG. 1. Degenerate spin wave modes in a collinear antiferromagnet with easy-axis anisotropy. 

Red and blue arrows represent the two sublattice spins Sa and 5'b in a unit cell, with z the easy-axis. (a) 
and (b): The two circularly polarized modes are characterized by left-handed and right-handed precessions 
around the easy-axis; they also have slightly different ratios between the cone angles of Sa and Sb- (c) 
and (b): The linearly polarized modes consist of different combinations of the circular modes. While 
Sa (Sb) is individually traveling counterclockwise (clockwise) on an elliptical orbit, the staggered field 
n = {Sa — Sb)/2S exhibits purely linear oscillation. 
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FIG. 2. A spin wave field-effect transistor and numerical simulation of its performance, (a) 

Schematics of the spin wave FET. Spin waves are generated by an oscillating magnetic field H{t) at one 
end of the device and then modulated by a gate voltage Vq. At the far end of the device, spin pumping 
into a heavy metal {e.g., Pt) induces an inverse spin Hall voltage which is measured by two voltmeters Vi 
and ¥ 2 - (b) Spin wave spectrum in the presence of the Dzyaloshinskii-Moriya interaction, (c) Numerical 
simulation of an amplitude shift keying on a AFM nanostrip (see Method). Spin wave signal of 1.4 THz 
is modulated by a 20GHz square wave from Vo. The relaxation time r « 3.5ps. (d) Parametric plot of the 
spatial pattern of the Neel order for selected sites in and around the gate during the interval t G [40ps, 50ps). 



FIG. 3. Spin wave state on the Bloch sphere. A spin wave state labeled by (0, (j)) corresponds to an 
elliptical precession of the Neel order illustrated in the right panel. 



















































